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ABSTRACT 

We have observed the Crab pulsar with the Deep Space Network (DSN) Goldstone 70 m 
antenna at 1664 MHz during three observing epochs for a total of 4 hours. Our data analysis has 
detected more than 2500 giant pulses, with flux densities ranging from 0.1 kJy to 150 kJy and 
pulse widths from 125 ns (limited by our bandwidth) to as long as 100 fj,s, with median power 
amplitudes and widths of 1 kJy and 2/is respectively. The most energetic pulses in our sample 
have energy fluxes of approximately 100 kJy-/is. We have used this large sample to investigate a 
number of giant-pulse emission properties in the Crab pulsar, including correlations among pulse 
flux density, width, energy flux, phase and time of arrival. We present a consistent accounting 
of the probability distributions and threshold cuts in order to reduce pulse- width biases. The 
excellent sensitivity obtained has allowed us to probe further into the population of giant pulses. 
We find that a significant portion, no less than 50%, of the overall pulsed energy flux at our 
observing frequency is emitted in the form of giant pulses. 

Subject headings: pulsars: general — pulsars: individual (Crab pulsar) 



1. Introduction 

While giant pulses (CPs) hav e been reli- 
ably detected from seven pulsars, (jKnighti 12006 
and references therein) their properties have 
only been well studied from two objects, the 
Crab pulsar, PSR B0531 -H 21 (ISallmen et al 



19991: ICordes et al.l l2004t iHankins et all l2003t 



Hankins &: Eilekl l2007l). and PSR B1937 -H21 



(jCognard et al.l 119961: ISoglasnov et al.ll2004D . In 

particular, the Crab pulsar has long been known 
as a giant pulse emitter. Its initial discovery 
(jStaelin fc Reifensteinll 19681 ) and a number of sub- 
sequent studies have reported remarkable proper- 
ties of giant pulses from the Crab. Giant pulses 



are a broadband phenomenon (e.g.. ISallmen et al 
19991) exhibiting the Galaxy's largest ob served 



brightness temperature (jCordes et al.ll2004[ ). and a 
subset of them are, in effect, super positions of ex- 
trem e ly narrow nanosecond pulses (IHankins et al 



2003: Hankins fc Eilek 200 



has been shown 
that energy flux emission from giant pulses exhibit 



opuise 
li. It 



^Jet Propulsion Laboratory, California Institute of 
Technology, 4800 Oak Grove Dr., Pasadena, CA 91109. 



a power law distribution jV(j!^ > Ep) ^ E^°' , wit h 
a in the rang e of 1.5-2.5 (IPoDOv fc S tapDerJ 2007l) . 
in contrast to normal radio pulse emission bein g 
Gaussian distributed (jHesse fc Wielebinskilll974[) . 

In this paper we analyze high time-resolution 
CP observations from the Crab pulsar conducted 
with Goldstone's Deep Space Network (DSN) 70 m 
radio telescope at 1.7 GHz (L-band). With our 
large data sample, we carried out a number of sta- 
tistical studies of CP properties. We present ra- 
dio observations and a description of the recording 
strategies and setup in Section 2. In section 3 we 
present the data reduction and analysis scheme. 
In section 4 we provide a description of the sta- 
tistical studies, followed by a discussion of the im- 
plications of the analysis. Finally, conclusions are 
given in section 5. 

2. Observations and Data Set 

As part of an effort to revitalize the L-band 
system on the DSN 70 m antenna in Goldstone, 
a number of pulsar observations were carried out 
over four epochs in early to mid 2008. The Crab 
pulsar was observed during two of these epochs 
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for a total observing time of almost four hours. 
For these observations, the front-end electronics 
down-converted the L-band RF signal to IF (via an 
intermediate S-band up-conversion) for recording, 
using a pair of VLBI Science Receivers (VSRs). 
The VSRs filter and sample the analog IF signal, 
then digitally form sub-channels for recording to 
disk. The data used for this study consist of 32 
MHz recorded bandpass in the range 1648-1680 
MHz, recorded as four adjacent 8 MHz channels 
with 2-bit I and 2-bit Q samples for 8050 continu- 
ous seconds, followed 3 min later by two 16 MHz 
channels with 1-bit I/Q sampling for 1830 continu- 
ous seconds. A number of hardware problems were 
encountered during the experiment's first hour, in- 
cluding antenna pointing and recorder sampling 
errors. The data corresponding to these problems 
were removed. The configuration of the channels 
is summarized in Table 1. 

The system temperature during Crab pulsar ob- 
servations may be dominated by the emission from 
the Crab Nebula, one of the brightest radio sources 
in the sky. This is certainly true in our case, 
where the nebula is not resolved by the antenna 
beam. The nebula is an extended source with a 
~ 5.5' diameter and a flux density paramet rized as 
Sn ^ 955i/-"-27 Jy (jBietenholz et al.lll997t ). where 
p is the observing frequency in GHz. At an ob- 
serving frequency of 1.7 GHz, the nebular flux 
density is Sn ^ 830 Jy. Since the DSN's 70 m 
antenna beam at L-band has a half-power width 
of 8.8', the nebula is not resolved so the sys- 
tem noise must correctly account for the nebular 
noise contribution. The nebular contribution Sn 
is combined with the contribution from the system 
temperature in absence of the Crab Nebula (S^y^) 
to obtain the total system temperature: 



Ssys — Sn 



S 



sys 



(1) 



On-off measurements of a standard calibrator 
3C48, and bright radio pulsar PSR B0329-h54 were 
carried out prior to observing the Crab. These 
observations yield a system temperature, S^yg, of 
35 K (10% error). Because the nominal gain G 
of the 70 m antenna is ~ 1.0 KJy~^, the mea- 
sured system temperature translates into a sys- 
tem equivalent flux density {Ssys — Tsys/G) of 
35 Jy. Adding the estimated nebular contribu- 
tion of ^ 830 Jy, yields a total system equiv- 



alent flux density Ssys = 865 Jy. This value 
agrees well with our estimates obtained from on- 
off measurements of the Crab, where we measured 
the system temperature while alternately point- 
ing at the Crab, and 1° away. We estimate the 
error for our overall flux density calibration scale 
to be less than ^ 20%. Our detection sensitiv- 
ity for single pulses is determined using the ra- 
diometer equation AS" = TjSsys/V Ai^At, where rj 
is the digitization loss factor, = 32 MHz is 
the recorded bandwidth and At = 0.1 /xs is the 
minimum sample time. With two-bit digitization 
we have 77 = 1.3 and we obtain a la detection 
threshold of ASmin 560 Jy. For a 7a detection 
threshold the minimum single-bin detectable pulse 
amplitude is expected to be 3.9 kJy. For compari- 
son. Table 2 lists the parameters of previous Crab 
GP studies with the current work. 

3. Data Reduction 

Our 4 hours of Crab pulsar data, consisting of 
540 Gigabytes, were recorded to disk and shipped 
to JPL for post-processing. The average and RMS 
RF voltage was computed for each second of data 
as a quick assessment of data quality. 

3.1. Coherent Dedispersion and Normal- 
ization 

The data were coherently dedispersed using 
the nominal dispersion measure (DM) valueQ of 
56.7671 pccm~^, followi ng the disp ersion removal 
technique developed bv. Hankins fc Rickett (1975|). 
Each IF channel was dedispersed separately by 
performing an FFT on a full second of complex 
(I and Q) samples, plus the fractional second 
of data required to fill the FFT arrays with a 
power of two number of samples. For the 2- 
bit, 8 MHz data, the approximately 1 Hz fre- 
quency bins were counter-rotated in phase to re- 
move the frequency-dependent dispersive delay, 
then inverse-transformed back to the time domain. 
The 1-bit, 16 MHz data were processed similarly 
except the final inverse FFT was performed on 
each half of the channel bandpass separately, ef- 
fectively splitting each 16 MHz channel into two 
8 MHz channels. In this way, all data could be 



^The DM value was obtained from Jodrell Bank Crab pulsar 
monthly ephemeris: 

|http://www.jb.man.ac.uk/~pulsar/crab.html[ 
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Table 1 
vsr channel configuration 



Epoch Number 


NRec^ 


Chan 1 


Chan 2 


Chan 3 


Chan 4 


BW (MUz)^ Nbits'^ 


Obs Time (min) 


1 


2 


1652 


1660 


1668 


1676 


8 2 


135 


2 


1 


1656 


1672 






16 1 


30 


3 


2 


1652 


1660 


1668 


1676 


8 2 


75 



^Number of recorders 
'^Channel Bandwidth 
'^Number of recorded bits per sample 



Table 2 

Crab pulsar giant pulse observations 



uiGHzY 


Epoch (MJD) 


T (hr) 


BW (MHz) 


AT ips) 


SEFD^ 




Telescope 


Reference 


0.1-0.2 


53635 




6 


1024-256 


1100 


31 


MWA-LFD 


(Bhat et al. 2007") 


0.430 


52304 


1.0 


12.5 


128 


1262 


11880 


Arecibo 


fCordes et al. 2004) 


0.6/1.4 


50224 




50 


1 




29 


vla 


fSallmen et al. 1999) 


0.812 


48433 


100 


20 


200-300 


13.5 


3 X lO'' 


GB 43m 


(Lundsren et al. 1995) 


1.18 


52277 


0.47 


100 


100 


309 


863 


Arecibo 


fCordes et al. 2004) 


1.197 


52944 


3.5 


20 


4 


215 


17869 


WSRT 


(Popov & Stappers 2007) 


1.475 


52277 


0.58 


100 


100 


291 


647 


Arecibo 


(Cordes et al. 2004) 


1.3-1.47 


53736 


3 


64 


0.5 


1100 


706 


ATCA 


(Bhat et al. 2008) 


1.7 


54618 


3 


32 


0.1 


865 


2500 


Goldstone 


This paper 


2.15 


52304-52306 


0.15 


100 


32 


79 


135 


Arecibo 


(Cordes et al. 2004) 


2.33 


52315 


0.15 


100 


32 


78 


92 


Arecibo 


(Cordes et al. 2Q04) 


2.85 


52306 


0.26 


100 


32 


74 


103 


Arecibo 


(Cordes et al. 2004) 


3.5 


52398-52412 


1.27 


100 


64 


41 


549 


Arecibo 


(Cordes et al. 2004) 


4.15 


52295-52337 


1.49 


100 


32 


20 


1663 


Arecibo 


(Cprdesjjt al. 2004) 


4.5-10.5 


53005-53736 




2200 


0.0004 




380 


Arecibo 


(Hankins k Eilck 2007) 


5.5 


52336-52411 


0.3 


100 


32 


20 


22 


Arecibo 


(Cordes et al, 2004) 


8.8 


52398-52414 


1.42 


100 


16 


22 


2249 


Arecibo 


(Cordes et al. 2004) 



^For comparison convenience this table is sorted by lowest observing frequency in ascending order 
''System equivalent flux density (Jy) 
'^Number of detected CPs 
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further processed in a similar manner. Dedispers- 
ing a full second of data kept delay smearing over 
each ~ 1 Hz frequency bin much less than a sam- 
ple, ensuring the dedispersed data maintained its 
full time resolution. 

After dedispersion, we form each sample's nor- 
malized power, Pi, averaged over all four frequency 
channels, where i is the sample number. This was 
done by independently normalizing each complex 
component (I and Q) from all four frequency chan- 
nels so that all eight quantities had zero mean and 
unit standard deviation. These eight quantities 
were squared and averaged to form our power time 
series, using: 



Q 



(2) 



k=l 



where k runs over the four frequency channels. As- 
suming each I and Q component is normalized and 
Gaussian distributed, in the no-signal limit Pi will 
have a reduced chi-squared distribution with eight 
degrees of freedom. 

3.2. Pulse Detection 

Our pulse-detection algorithm begins with the 
8 MHz Pi time series whose noise component is 
nominally modeled by the reduced chi-squared dis- 
tribution with 8 degrees of freedom. In general, 
with V degrees of freedom, the reduced prob- 
ability density function, parameterized by x, is 
given by: 



pdf(a;) = 



i>{i'x) 



e 2 



2'^/2r(j./2) 



(3) 



where the F function for integer and half-integer 
arguments is defined as: 



^(„) = ,»-l)^ r(i„) = (4) 

The red data points in Figure [1^ show the mea- 
sured Pi distribution for 30 minutes of data. The 
data samples were taken during the Crab pulsar 
off-pulse regions of the pulse phase. The black 
solid curve in [T^ is the corresponding theoretical 
expectation given by Equation |3] with i/ = 8. Ex- 
cellent agreement is seen in over eight decades. 



In the time domain, our matched filter algo- 
rithm averages iV temporally consecutive Pi power 
samples, and compares this with a threshold pa- 
rameter designed to strongly reject noise while 
passing GPs with high efficiency. The filter used 
in this analysis begins with the Pi power series 
corresponding to N=l. The N=2 time series is 
formed by averaging each consecutive pair (with- 
out overlap) of the N=l series, resulting in half the 
number of power samples. This doubling proce- 
dure is repeated up to the creation of the 7V=2048 
power series which corresponds to 256 lis wide 
time-averaged bins. For each power series, sum- 
ming over the N samples results in: 



1 



N-l 



m— 



(5) 



where i is in the range (0, 2N, ...). Because each 
Pi averages 8 Gaussian squares, Pf^ will have its 
noise component distributed according to Equa- 
tion [3] with v = 8N, assuming the data samples 
Pi are truly independent^ Figure [1^ shows the 
off-pulse power, Pf^ , for = 1 (red points ly — 8 
in Equation |3|) , N = 2 (green curve, v = 16) and 
N = A (blue curve, v = 32). The black curves 
give the theoretical distributions from Equation |3] 
in each case. 

In order to prevent cut-induced biases in the 
giant-pulse width distribution, proper threshold 
cuts should result in the same number of back- 
ground noise events per second, on average, re- 
gardless of A^. We define R to be the rate of noise 
fluctuations expected to pass a threshold cut ^{N), 
in units of events per second. This is equivalent 
to picking a constant event confidence level inde- 
pendent of pulse width. The noise rate R, also the 
false alarm rate, is equal to the probability that an 
averaged power computed with A^ samples exceeds 
^, times the number of trials in 1 second: 



(6) 



where 5(^) is 1 minus the cumulative probabil- 
ity distribution corresponding to Equation [3J also 



^As A'^ increases, a small sample-to-sample correlation was 
seen. This loss of true data independence was handled by 
making a small correction to the ideal v. For example, the 
Ns = 4 case, ideally i/ = 32 but u=31 yields a better fit to 
the data. 
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Fig. 1. — a) The background pulse xl distribution for of 8, 16 and 32; the red, green and blue data points 
respectively, b) The background pulses plus giant pulse xt- The corresponding theoretical xt distributions 
are the solid black curves, c) The noise rate survival distribution for v of 8, 16 and 32; the red, green and 
blue curves respectively. The black horizontal dashed line, at a noise rate of 1 Hz, intersects the three Xu 
distributions at the corresponding power cut points. Three black vertical dashed power cut lines have been 
extended through all three figures. 
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known as a survival function, and fs is the num- 
ber of samples per second. Figure [TJ; shows the 
the noise rate as a function of P for the = 1 
(red points), N — 2 (green curve, = 16) and 
iV = 4 (blue curve, v — 32) cases. 

In the limit of large ^, reasonable given our de- 
sire to strongly cut noise, this becomes: 



R = 



4Af-lg-4Af5 



N 



(AN - 1)! 



(7) 



The threshold values are computed from this equa- 
tion iteratively using 



4^[(4iV-l)log(0 

where 



log(i?)-t-C] (8) 



C={m- 1) log(47V) + log(^) - log((47V - 1)!) 

(9) 

For each averaged power P.^ in every power se- 
ries. Equations [7] and [8] are used to compute the 
effective noise rate R corresponding to the mea- 
sured Pj^. This effective noise rate estimates the 
likelihood that a given averaged power value is a 
noise fluctuation. The dotted horizontal line in 
Figure[TJ; is the 1 Hz rate cut, i?-cut=1.0 Hz, and 
the corresponding power cuts ^{N = 1), S^{N ~ 2) 
and £,{N = 4) are shown as the three vertical dot- 
ted hues in Figure [1^ and Figure [T}). As described 
below, we reduce background noise by cutting on 
R, using the £,{N), in order to guarantee that the 
background rate is independent of N. 

This consistent accounting of the probability 
distributions and equivalent threshold cuts to keep 
background rates independent of pulse width en- 
sures balanced detection efRciencies even for ex- 
tremely large- width giant pulses. 

Figure [It plots the xl distributions for all the 
data, so as to include the Crab's giant pulses: the 
long tail of events with large xt is clearly seen. 
With each detection, both A'^ and the probability 
of such a pulse being a noise fluctuation, expressed 
as R in Equation [3 are stored. 

This pulse-detection algorithm typically detects 
each pulse candidate a large number of times, so 
a method is needed to transform these multiple 
detections into individual pulse candidates, and 
to estimate their properties such as time, total 
power, and width. The algorithm used in this 



analysis was to flag any sample that participates in 
a detection, regardless of N, then call any contigu- 
ous interval of flagged samples a pulse candidate. 
A number of parameters were computed for each 
pulse candidate. The peak power, the pulse inte- 
grated power, the time at peak power, the power- 
weighted mean pulse time, and the power variance 
about the mean time, were all computed from the 
8 MHz channel- summed powers. From the mul- 
tiple detections associated with each pulse can- 
didate the lowest-probability R, corresponding to 
the highest significance, was chosen as the giant 
pulse. The associated N, power, and time were 
then saved. 

Figures [5^ and show scatterplots of pulse 
phase and time for all giant pulses found in our 
data set, for both a relatively soft (i?=0.02 Hz) 
and hard cut (i?=0.0006 Hz), respectivelsH. The 
main pulse, with phase near 0.3, and the intcr- 
pulse, near 0.7, are clearly evident. The loose R- 
cut (0.02 Hz) implies more white-noise contami- 
nation due to the corresponding lower xt thresh- 
olds. Similarly, the tight R-cut (0.0006 Hz) re- 
duces noise contamination on higher x^, which is 
equivalent to a higher power signal-to-noise (SNR) 
threshold. The projections of these scatterplots 
onto the time axis, the pulse phase histograms, 
are presented in Figures [2j: and [2ji for the hard 
and soft cut, respectively. In the loose-cut case, 
over 1600 main giant pulses (GPs) are seen on a 
fiat background of 100 noise pulses. With tight 
cuts, over 1200 main pulse GPs are found with 
no background, and a clear peak of GPs is seen 
at the interpulse phase region. By defining an on- 
pulse phase range of (0.32, 0.35) for the main pulse 
phase, and (0.72, 0.75) for the interpulse phase, we 
can study the count rate and SNR of the GPs as 
a function of i?-cut. 

If St is defined as the number of counts within 
the on-pulse region and Sn is the number of back- 
ground counts estimated from the off-pulse phase 
region, the background-subtracted signal count is 
S = St — Sn- Assuming counting statistics, the 
SNR is then: 



SNR^{ST~SN)/iST + SNy^^ (10) 



^The significance of these cut values becomes clear later in 
the paper. 
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Fig. 2. — Time phase scatterplots of the detected pulses for the loose and tight data sets. The main and 
interpulse giant pulses are seen as straight lines at phases of 0.337 and 0.737 respectively. The bottom 
histograms shows the corresponding phase projection. The improvement in signal to noise is clearly seen 
with the tight cuts. 
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Figure [3] shows the total number of giant pulses 
as a function of i?-cut, for both the main and in- 
terpulse regions. For the loosest cut processed 
(i?=0.6), over 2500 main phase GPs and 200 inter 
phase GPs are observed. 
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Fig. 3. — The number of background-subtracted 
Giant Pulses, for both the Main pulse and Inter- 
pulse, as a function of the noise cut, i?-cut. As the 
i?-cut increases (power amplitude cut decreases) 
the total number of giant pulses increase. 

Careful examination of the pulse-candidates in 
Figure [2^ shows density variations during the ex- 
periment's first and last hours. Our detection 
thresholds were adjusted to allow an approxi- 
mately constant noise rate as a function of pulse 
width, however Figure [5^ shows that the noise 
rate is not constant in time. At the experiment's 
start, there were a number of problems with the 
front-end sampler voltage thresholds, which re- 
duced the data sensitivity. The single-bit data 
beginning at time 8230 s also shows a low pulse- 
candidate density. Even though dispersion mit- 
igates the 1-bit sampling loss to some extent, as 
explained above, these data are somewhat less sen- 
sitive than the 2-bit data. Efforts to reconcile the 
different sampling effects and resulting sensitivi- 
ties, for example by scaling the sample values or 
modifying the detection thresholds, significantly 
increases the complexity of the analysis. Thus, 
from this point forward the analysis is restricted 
to only those data taken from time 2400 to 8050 s. 



a total of 5650 seconds of data, as these data have 
the highest sensitivity, and show an approximately 
constant sensitivity. 
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Fig. 4. — The giant pulse signal-to-noise ratio as a 
function of the noise cut, i?-cut, for both the main 
pulse and interpulse GPs. 

Figure |4] shows the SNR as a function of i?-cut 
for both the main and interpulse pulses. It is seen 
that the maximum SNR for the main pulses are 
at i? = 0.02 and i?=0.0006 respectively, lower- 
ing the i?-cut (increasing the total cut) fewer 
GPs with higher SNR values are detected. In this 
paper we will refer to three cuts, a loose cut at 
i?=0.02, a tight cut at i?=0.0006, and the loos- 
est cut at i?=0.6 resulting in the largest sample 
of Main Pulse (MP) GPs found in this analysis 
(~2500). Yet even with the tight cut over 1200 
MP GPs and over 100 Inter-Pulse (IP) GP are 
detected both without any significant background 
contamination. 

4. Various Statistics of Detected Giant 
Pulses 

The tremendous amount of energy radiated by 
giant pulses is observed with a wide variety of 
pulse morphologies. One of the goals of giant- 
pulse research is to determine, as closely as pos- 
sible, the originating pulse shape, its amplitude 
and width. A wide variety of factors influence 
the shape of the pulse, such as dispersion, in- 
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strumental smearing, interstellar scintillation and 
scattering due to turbulent media. These effects 
all tend to broaden a pulse, decreasing its true 
peak amplitude and increasing its width. In our 
data set, the amplitude is simply defined as the 
peak power, or rather the peak flux density, found 
in the optimum smoothed data set. Past attempts 
to characterize pulse widths for a large sample of 
GPs have been limited either to coarse sampling 
size or l arge dispersive smearing and small sam- 
ple sizes j Lundgren et al.|[l995 ; Popov fc Stappers 
[2007; B hat et al.ll2008l ). Our L-band observations, 
employing coherent dedispersion on large band- 
width baseband channels minimize the dispersive 
smearing, obtaining a 125 /iS time resolution, con- 
tain over 1200 GPs, which should allow more ro- 
bust width measurements. 

4.1. Pulse Amplitudes and Widths 

Visual examination of many pulses, some of 
which are shown in FiguresE^-f, reveals that varia- 
tions in pulse morphology represent the dominant 
systematic error in pulse amplitude and width de- 
terminations. The six GP events shown are plot- 
ted using their optimum smoothing, with widths 
ranging from the narrowest (smoothing widths of 
iV=4) to the widest (A^=512), and with ampli- 
tudes ranging from 1 kjy to above 100 kJy. Al- 
though a large subset of GPs can be well fit to 
a Gaussian shape, an equally large sample of GPs 
show non-Gaussian shapes, as presented in Figures 
[St-f. Given a significant fraction of GPs having 
non-Gaussian shapes, fitting GPs with a Gaussian 
model yields poor chi-squared fits, and thus poor 
amplitude and width estimates. 

An alternative width definition is simply the 
optimum smoothing width Wg = NAt {At is the 
intrinsic sample resolution) obtained from the gi- 
ant pulse detection algorithm. The amplitude and 
timing can then be defined using the bin with the 
greatest power amplitude. With these alternative 
definitions of power amplitude and width, pulse 
morphology can be seen in Figure [6^, where the 
scatter plot shows GP peak flux density versus 
width, for all the GPs passing the tight-cut. A 
strong correlation between the observed peak flux 
density and smoothing width is observed, with 
higher peak flux density seen at lower smoothing 
widths. The diagonal dotted lines in Figure [6^ 
show iso-energy flux contours of 1 kJy-/is, 10 kJy- 



fis and 100 kJy-/Lts. A wide variation of shapes 
can be seen for the same observed energy flux. 
For example, a large pulse with E = 10 kJy-/is 
can be very sharp, less than a /is in width with 
over 10 kJy peak flux density, or very broad, with 
widths over 10 /iS but flux density less than 1 kJy. 
Although great care was taken to ensure that no 
pulse width selection bias occurred, a comparison 
of the GP data with the iso-energy flux contours 
suggests that giant pulses with total energy fluxes 
of ~1 kJy-/j,s cannot be extracted above smooth- 
ing widths of 2 lis. 

Figures [6}d and Eh show the corresponding am- 
plitude and width projections. The peak flux den- 
sity distribution in Figure Et) indicates the median 
pulse flux density is ~1 kJy with a minimum ob- 
served peak flux density of ~0.1 kJy and a maxi- 
mum of ~100 kJy. Above a peak flux density of 2 
kJy, the histogram shows a power-law dependence 
out to 100 kJy. A power-law fit yields an expo- 
nen t of a = -2.2 +/- 0.1, very close to that found 
by (jBhat et al. I l2008h . where at slightly lower fre- 
quencies of 1300 MHz and 1460 MHz, values of 
-2.3 and -2.2 were measured respectively. 

The distribution of the smoothing width, Ws , is 
shown in FigureEfc- The median Wg is 2 /is and the 
mean Ws is ~6 /iS with a minimum of 0.2 /xs with 
a long tail extending out to 100 /iS. This width 
range compared to previous Crab GP datasets 
( Bhat et al. 20081) is extended to pulses with 10 
times wider and two times narrower widths. 

Although the smoothing width, Wg is a conve- 
nient definition of pulse width, it does not always 
appear to be optimum. Sometimes, as shown in 
Figure \5jp, a wide giant pulse with A^=512 has 
almost all its energy within the one smoothed 
bin. Yet as Figures Et-f illustrate, a large frac- 
tion the time signifiant energy is spread over many 
bins at the optimum smo othing width. Following 
others ( Bhat et al. I l2008h we define an "effective 
pulse width," Wgtd, which is the standard devia- 
tion (std) for all "acceptable" bins in a giant pulse, 
where acceptable bins are defined with a simple 
algorithm. Raw data with the highest time reso- 
lution is used, and data within ~3A'^ of the max- 
imum power amplitude sample is selected. Only 
bins with power amplitudes greater than 2 sigma 
above the background power are included in the 
first estimation of the Wstd, and the background 
level. The final calculation sums over all bins 
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11 



within '--^SWstd of the peak, again with bins having 
power amphtudes greater than 1 sigma above the 
background. This algorithm allows us to estimate 
the pulse standard deviation and higher order mo- 
ments such as skewness and kurtosis. This proce- 
dure gives an estimate of the pulse width (Wstd) 
that, as expected, gives excellent agreement with 
the fitted Gaussian sigma and gives a reasonable 
estimation of the pulse width in the non-Gaussian 
cases. A comparison with pulse widths extracted 
by Bhat et. al., at 1470 MHz, can be seen in Fig- 
ure [71 Due to the higher time resolution in this 
study, pulses with widths a factor of 2 narrower are 
detected. Perhaps the more significant difference 
with the earlier data of Bhat et. al. is the increase 
in the sensitivity at very large pulse widths. Bhat 
et. al. observed no pulses with widths greater 
than 10 ^s, whereas the present data set has a 
significant number of pulses extending out to over 
100 /xs with a median Wstd of 2.35 /iS and a mean 
of 10.45 /zs. We have parameterized this distribu- 
tion with a good power law fit, N{Wstd) W^Tj^'^, 
that extends from ^ 1 fis to almost 100 fis. 
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Fig. 7. — The giant pulse width histogram for the 
tight data set. The fitted power law is indicated 
by the solid line. The dotted histogram is data of 
(iBhat et al.ll2007t) 



4.2. Pulse Energies 

The measurement of giant pulse parameters 
such as power amplitude, width, and energy flux 



are critical in constraining models of GP emission. 
The total pulse energy flux has the advantage, as 
compared to pulse width, of being a robust pulse 
estimation parameter, as it is relatively insensitive 
to pulse shape. Also, since the total pulse energy 
flux is a simple power sum over the pulse duration, 
the random background noise will tend to average 
down to zero. The GP energy flux distribution, 
for the tight cut data set, is shown as black circles 
in Figure [51 Energy fluxes up to ^100 kJy-/xs are 
found, consistent with earlier studies ( Bhat et al 



The average observed GP energy flux for 
all the main giant pulses was ~7 kJy-^s with a 
median energy flux of ^^5 kJy-/xs. Above 5 kJy-/is 
the data are well fit by a power law, the solid line 
in Figure [H shows fit to the data with q;=-2.57. A 
very sharp energy flux threshold is observed at ^3 
kJy-/is. The black diamonds in Figure [51 show the 
main GP energy flux distribution for the loosest 
cut examined. An additional ~1400 main GPs are 
found but all at lower energy flux as would be ex- 
pected. The additional main GPs seem to follow 
the same power law as the data with tight cuts 
but the energy threshold is reduced to ~1.5 kJy- 
fis. No evidence is seen for a roll over or softening 
of the intrinsic power law. The turnover in Figure 
[His most likely due to the applied threshold cut. 

Following a previous analysis (Knight 20061 ) . 
the cumulative energy flux distribution, deflned 
as the probability per second of a pulse having an 
energy flux E greater than a value Eq, is used to 
compare the occurrence frequency of giant pulses 
between different experiments: 



P{E > Eq) = kE°' 



(11) 



Figure [9l shows the observed cumulative energy 
flux distribution for the tight data set, solid black 
circles, along with the data from Bhat et. al., open 
squares. The solid line is a power-law flt to the 
data in the energy range 5 to 100 kJy-/zs. Above 
the energy flux threshold, the data are well flt by a 
power law and no evidence is found for high-energy 
flux cut-offs. The best-flt power law exponent is 
found to be a=-1.9 +/- 0.05 The agreement with 
Bhat et. al. is very good in the energy flux region 
of 10-20 kJy-/xs but significant differences in slope 
and absolute magnitude are seen at higher and 
lower pulse energy fluxes. This disagreement can 
be understood by a combination of differences in 



12 



1000 



100 



- •- - Tight Cuts 

Power law fit a--2.57 
- Loosest Cuts 



cn 
c 

O 

O 




10 - 



Although earher work by Bhat et al. (jBhat et al 



1 10 100 1000 

Energy (kJy-^s) 

Fig. 8. — The giant pulse energy flux distribution 
with a power law fit of a = -2.57 for energy fluxes 
between 5 and 50 kJy-/is. 



antenna sensitivity and pulse algorithm efficiency. 
As discussed earlier, for any fixed pulse energy a 
wide range of pulse amplitudes and widths are ob- 
served. A loss of pulse detection efficiency at low 
and high widths may alter the true nature of the 
cumulative energy flux distribution. 
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2008f) only co vered widths from 1 to 10 ^s, others 
(jPopov fc Stappers 2007) found a dependence on 
pulse width, such that at narrow widths (^4 /is), 
a=-1.7 whereas for wider pulses (65 /is), a was 
found to be =-3.2. To examine this dependence 
our tight data set was divided up into two sub- 
sets: a narrow {Ws less than 5 /is) and a wide 
data set {Ws greater than 5 ^s). In both cases, 
the slope of the cumulative energy flux distribu- 
tion was estimated. The narrow data set yields 
a slope values of -1.7 to -2.0, dependent on the 
energy flux range fit, which is very consistent with 
Bhat et al. While the wide data set yields a steep 
slope of a~-2.4: to -3, consistent with the work of 
Popov and Stappers. 

It is interesting to compare the total energy flux 
emitted by the Crab pulsar as giant pulses with 
the overall pulse emission at this frequency. Us- 
ing the background- free data set (tight cuts), we 
summed the energy flux of each candidate GP, ob- 
taining an average pulse energy flux of 0.134 kJy- 
Ijs per rotation over the span of the observation, 
which amounted to ^ 168,000 rotations. On the 
other hand, the mean pulse energy flux from the 
pulse profile containing all rotations is 0.26 kJy- 
fj,s. This value is in excellent agreement with the 
quoted value in the ATNF catalog after extrap- 
olating the mean flux density from 1400 MHz to 
1665 MHz using their spectral index of -3.1. This 
suggests that a significant fraction, at least ~50%, 
of the pulsar emission energy may be emitted as 
GPs. Another estimate for pulse energy fiux ob- 
tained at a frequency of 1664 MHz, nearly iden- 
tical to our observing frequency, suggests a mean 
energy flux of 0.125 kJy-^s for the main pulse and 
0.025 kJy-/xs for t he interpulse at this frequency 
( Mancheste 3l97lfl . The total energy flux amount- 
ing to 0.15 kjy-yus is very close to the entire GP 
energy flux emission, further suggesting that a sig- 
niflcantly large portion of the energy flux may be 
emitted as GPs. 

As a further check we compared the pulse pro- 
file including all rotations with a second profile, 
where rotations containing an identified GP is re- 
moved from the profile. Figure ITOl shows the two 
pulse profiles. The second profile (dashed curve) 



■^ATNF pulsar catalog web address: 
|http : / / www, atnf.csiro.au/ research / pulsar / psrcat | 
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is obtained using a loose cut to identify GP can- 
didates. The figure shows that GPs make up 54% 
of the overall pulse energy flux at our observing 
frequency. Also, evident from the figure is how 
identical the two profiles are in terms of the pulse 
width. At the moment we do not have sufficient 
statistics to resolve and quantify the width of each 
profile and instead rely on a rough visual estimate. 
This is further indication that a similar significant 
reduction in energy flux is obtained assuming that 
the GP and main pulse widths are similar. 
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Fig. 10. — Crab pulse profiles obtained (solid 
curve) using the entire sample discussed in this 
analysis. The dashed curve shows the profile with- 
out rotation cycles that include a GP candidate. 

A large fraction of the emission energy flux at 
this frequency is already accounted for with GPs 
identified at our current sensitivity. If the power 
law continues to lower energies, further improve- 
ment in sensitivity will allow us to reach the in- 
trinsic turnaround in the distribution necessary to 
keep the total GP emission energy below the total 
pulsed energy. In this case, there would be only 
one fundamental process responsible for the total 
pulsed emission in the Crab pulsar. 

4.3. Pulse Asymmetry and Shape 

Skewness can be used as a measure of GP asym- 
metry, and is defined here to be: 

s = s]Ty,(t, (12) 



where the weighted sum is over a narrow region 
of the pulse with respect to its peak position tp. 
The weight, Wi, is the power amplitude level in 
the corresponding time bin, cr is the standard de- 
viation (yVstd) and the pulse algorithm used is the 
same as that used for obtaining the pulse standard 
deviation. For a Gaussian shaped pulse, skewness 
will have a value near zero. Pulses with long tails 
following the peak will have a positive skewness 
and those with tails before the peak will have neg- 
ative values. 

Figure [TT] shows the skewness distribution for 
both the giant pulses and background pulses, solid 
and dashed histogram respectively. Both distribu- 
tions show average skewness that is very close to 
zero with very similar spread in the distribution. 
No statistically strong skewness dependence was 
found as a function of width, amplitude or energy 
and no differences were found between the inter- 
pulse and main giant pulses. 



350 




Skewness 



Fig. 11. — Skewness histograms for both the giant- 
pulses and background pulses for the tight data 
set. 

The next higher statistical moment, kurtosis, 
can also yield useful information on pulse shape, 
and is defined as: 

k = ^W^{U-tpf/{^W^)cJ^-i. (13) 

Kurtosis is zero for a Gaussian-shaped pulse, neg- 
ative for pulses that are flatter than Gaussian, 
and positive for pulses more steeply peaked than 



14 



Gaussian. To account for algorithmic dependences 
and background noise, we have calculated skew- 
ness and kurtosis for both giant pulses and out-of- 
phase background pulse-candidates. 

Our kurtosis distribution is shown in Figure 
[T2I for both GPs and background events. The 
background pulses are peaked near zero, consistent 
with the expectation for gaussian shaped pulses. 
A positive kurtosis tail extending to about a value 
of 5 is observed. This most likely reflects the bias 
of the decimation portion of the matched filter de- 
tection algorithm, where smoothing tends to in- 
crease kurtosis. The Crab's giant pulses are seen 
to be peaked at a kurtosis value of 0.6 with a 
very long tail extending to a value of 50, imply- 
ing they are more sharply peaked than pulses from 
the white-noise background. The reason for this 
can be seen by examining the dependence of mean 
kurtosis on the total GP energy flux shown in Fig- 
ure [ini This plot shows GP averaged kurtosis in- 
creases linearly with the total energy flux. Al- 
though the background pulse candidates do not 
extend to high energies, it is clear their behavior 
is quite different from the GPs upward trend. 
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Fig. 12. — Histograms of pulse kurtosis for both 
giant pulses and background pulses for the tight 
data set. 



4.4. Pulse Time of Arrival 

The Crab pulsar has a well defined period of 
^33 ms and the phase of the main and interpulse 
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Fig. 13. — The dependence of kurtosis mean with 
the total pulse energy flux for the giant pulses in 
the tight data set. 

giant pulses are well regulated within the pulsar 
period. The phase of both the interpulse and main 
pulses are seen in the projections of the scatterplot 
in Figure [2j The pulse phase residual, or Time of 
Arrival (TOA) residual, is the TOA of the GP 
peak with respect the start of the pulsar model 
cycle. For the Crab, at L-band frequency, both the 
main and interpulse GPs are found to have TOA 
residuals that fall within 1% of a cycle 330 /is). 
Figure [14] shows the histogram for the main and 
interpulse GPs, where the mean TOA has been set 
to zero for both. To remove possible confusion in 
the interpulse data set due to background pulse- 
candidates, the tight-cut data set has been used. 
The main pulses are found to have rms of 0.0028 
± 0.0001 cycles or 1° of phase corresponding to ~ 
90 /is, while the interpulse GPs where found to be 
wider with an rms of 0.0042 ± 0.0003 or 1.5° of 
phase corresponding to ^ 140 /iS. 

Since it is expected that the intrinsic pulse 
phase jitter and smearing will be independent 
of absolute phase location, the 50% increase in 
TOA width is significant and may give an im- 
portant constraint in astrophysical pulsar mod- 
els. Two earlier studies, both with smaller gi- 
ant p ulse samples, ( Bhat et "al]l2008 : Cordes et al 



2004l) found evidence for stronger pulses to have 
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Fig. 14. — Histograms of the time of arrival for 
the Main and Inter-pulse Giant pulses for the tight 
data set. 
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Fig. 15. — A scatterplot of the phase and total 
pulse energy flux. 



narrower phase windows. Where strength was de- 
fined as the peak pulse amplitude. Due to the 
strong correlation of peak flux density and width, 
as shown in Figure 7a, it is reasonable to believe 
that the narrow widths of the larger peak pulses 
may indeed result in better phase resolution. To 
examine this possibility the tight data set was di- 
vided into two subsets, a large and small peak flux 
density data set, those with peak amplitude larger 
than 5 kJy and those smaller than 1 kJy. The 
widths of the TOA residuals for the two n GP sam- 
ples were then estimated to be cr= 0.0027 ± 0.0002 
and cr= 0.0027 ± 0.0001 for the large and small 
peak samples respectively. No significant differ- 
ence is found, as the sample TOA residuals are 
seen to be statistically consistent to each other. 

A plot of the main pulse energy fiux versus 
phase is shown in Figures [TSis scatterplot. No cor- 
relation between phase and energy fiux is evident. 
To examine this more quantitatively, the data were 
divided into three total-energy fiux classes: high 
energy fiux (E > 200 kJy-/is), medium energy flux 
(100 kJy-^s < E < 200 kJy-^s), and low energy 
flux (E<100 kJy-^s). The table below summa- 
rizes the TOA residual statistics of each data set. 
No tendency is seen for larger energetic pulses to 
originate in narrower phase windows. 

The joint statistics of GP total energy flux with 
the interarrival time (lAT) of the next GP is shown 



in Figure [T6l The lAT is simply the measured time 
between between one GP and the next. The in- 
terarrival time for both the Main pulse (MP) and 
Interpulse (IP) giant pulses are plotted, where the 
lAT is converted to the number of crab pulsar ro- 
tations (~ 33 ms). The discrete nature of the lAT 
for pulses below 10 rotations is simply a reflec- 
tion of the phase cut (0.30-0.35) for the MPs and 
(0.72-0.75) for the IPs. For the sensitivity achieved 
with our observations, with tight cuts, the average 
number of Crab rotations between giant pulses was 
found to be 127 with the largest gap between gi- 
ant pulses at ~ 1000 rotations. With the loosest 
cut used, the averaged number of Crab rotations 
between GPs are 67. No observable dependence 
between the total GP energy flux and the arrival 
time of the next pulse was seen. 

The smallest I AT of ~ 0.4 cycle was observed 
twice, these cases being a main pulse immediately 
preceded by an interpulse (I-M) . No cases of inter- 
pulse immediately following the main pulse were 
found in the tight data set (M-I), corresponding to 
an lAT of ~ 0.6 cycle. To examine the indepen- 
dence between main and interpulses, we examined 
the consistency of the data with the assumed Pois- 
son nature of the giant pulse process. A Poisson 
process describes events which occur continuously 
and independently, and for this study would imply 
the interarrival times are exponentially distributed 
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Table 3 

Giant Pulse TOA residual statistics. 



E ( kJy-Ms) 


Sum 


Mean 


St.Dev. 


E < 100 


564 


0.3368 +/- 0.0001 


0.0028 


100 < E < 200 


571 


0.3373 +/- 0.0001 


0.0028 


E > 200 


106 


0.3365 +/- 0.0003 
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Fig. 16. — A scatterplot of the Giant-Pulse energy 
and and Interarrival times for the Giant-Pulses in 
the tight data set. 



with parameter A (mean = 1/A). 

P{t) = Ae" 



(14) 



Figure [17] shows the distribution of MPs inter- 
arrival times along with the best fit to a Poisson 
model. The data appears to be consistent with 
a Poisson process, the best weighted fit yielding 
A=0.223 +/- 0.007 Hz with a x^=19 for v=2?,. 
This A value agrees well the inverse of the mean 
interarrival times, 0.234 +/- .007. Within statis- 
tics, it appears that the arrival times of Crab's gi- 
ant pulses are Poisson distributed, implying they 
are memory less, and any given time interval is in- 
dependent of what occurs before or after. 

In addition, both the interpulses and the com- 
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Fig. 17. — The interarrival time of main pulses 
(MPs) for the tight data set. The solid curve is an 
exponential fit to the data. 



bined GP data sets show interarrival time distribu- 
tions that are also consistent with a Poisson pro- 
cess. We conclude that interpulse GPs do not ap- 
pear to be correlated with main pulses to our level 
of examination. 

5. Summary and Conclusions 

We have examined the energy flux distribution, 
timing, and statistical properties of giant radio 
pulses from the Crab pulsar at 1700 MHz using 
dsn's 70 m antenna at Goldstone, achieving a 
time resolution of 125 ns. Our pulse detection was 
based strictly on the underlying Chi-square statis- 
tics of the data, where we attempted to keep the 
noise rate and event confidence levels constant and 
independent of pulse width. Our consistent ac- 
counting of the probability distributions and non- 
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Gaussian effects, in order to keep background rates 
independent of pulse width, is an improvement 
over previous analyses. 

The statistical analysis of giant pulse popula- 
tion was carried out using three significance cuts: 
a loose cut (R-cut = 0.0200), where we obtained 
1879 pulses, and a tight cut (R-cut = 0.0006) 
yielding 1314 pulses and the loosest cut studied 
(R-cut — 0.6000) which yielded over 2500 main 
GPs. With a large number of giant pulses we 
were able to study various statistical properties 
of GPs. We have confirmed the power-law na- 
ture of the peak pulse flux density distribution and 
have obtained a power-law slope of —2.2, consis- 
tent with (Bhat et al..,2008 ) results. We observed 
giant pulses with widths up to 100 ^s, 10 times 
larger than the maximum found by Bhat et. al. 
Further, the distribution of pulse cumulative en- 
ergy fluxes also follows a power-law distribution 
with a slope of -2.57 , consistent with Bhat et 
al. (jBhat et al.l|2008) . We also looked for cor- 
relations between the time of arrival of individ- 
ual giant pulses and concluded the time of arrival 
of each event is consistent with Poisson statistics. 
Our pulse detection sensitivity, high-time resolu- 
tion, and consistent statistical treatment of the 
data has allowed us to examine a weaker popula- 
tion of giant pulses, narrowing the gap with "nor- 
mal" pulses from the Crab pulsar. 

Our comparison of pulsed energy flux emission 
from GPs with the normal pulsed emission has 
forced us to claim that a significant portion, per- 
haps as large as ^-^90% but at least 50%, of the 
emission at this frequency is in the form of GPs. 
With further improvement in sensitivity and ob- 
serving time, it should be possible to either con- 
firm the turnaround seen in the energy distribu- 
tion or push it only slightly lower before reaching 
the intrinsic turnaround in the distribution. Such 
a study could be carried out by a phased array 
instrument such as the VLA, where the narrower 
beam will make it possible to lower the nebular 
contribution to the overall system temperature, 
providing better sensitivity to further probe the 
weaker population of GPs as well as the remain- 
ing normal pulses in the system. 
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